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Abstract 

We consider solutions to the complex Trkalian equation, V x c = c, where c is a 3 component 
vector function with each component in the complex field, and may be expressed in the 
form c = e* 9 V-F, with g real and F complex. We find, there are precisely two classes 
of solutions; one where g is a Cartesian variable and one where g is the spherical radial 
coordinate. We consider these flows to be the simplest of all exact 3-d solutions to the 
Euler's equation for the ideal incompressible fluid. Pictures are presented. 

The novel approach we use in solving for these classes of solutions to these 3-dimensional 
vector pdes involves differential geometric techniques: one may employ the method to gen- 
erate solutions to other classes of vector pdes. 
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1 Introduction 



The Beltrami equation 

vx(Vxv)=0 (1.1) 

has received enormous attention in recent years. Reviews§~! note the prominent role that 
Beltrami fields play in the theory of exact, closed form solutions to the Euler's and Navier 
Stokes equations and their relations to the electromagnetic wave equations. Moreover, 
Beltrami fields are related to minimum energy plasma fields and have therefore garnered 
much attention from the MHD community.@~0 In this paper, we would like to undertake 
the first steps in a systematic classification of these vector fields, using a novel differential 
geometric technique which we discss below. 

Any 3 component vector function, v, in 7Z 3 may be written as the real part of a three 
component complex vector field, of the following form: 

c = e i9 VF. (1.2) 

where g is a real function and F is a function from 7Z 3 to C 3 . Tr kalian fields! are real vector 
fields which solve 

Vxv = v, (1.3) 
and therefore satisfy 1.1. We solve the complex version of the Trkalian equation: 

V x c = c, (1.4) 

where c is also a vector of the form of Eq. (1.2). Notice that both the real and imaginary 
parts of c must solve the Trkalian equation, which is a linear equation. From solutions to the 
Trkalian equation, one straightforwardly may derive (non-zero helicity) three dimensional 
(in a topological sense: see ||) time- independent solutions to Euler's equations for the in- 
compressible ideal fluid0 and time dependent solutions to the incompressible Navier Stokes 



equations.0 Therefore we consider the equations (1.2) and (1.4) to be the simplest possible 
three-dimensional hydrodynamical equation. Although we do not discuss the stability of 
these solutions, a discussion of boundary conditions is given in section 4. 

We completely categorize all solutions to the simulataneous set of equations given by 
(1.2) and (1.4). The two classes of solutions to this problem are A.) where the function g 
in 1.4 is identical to a Cartesian variable (for example: z) and the function F is a complex 
analytical function of the other two Cartesian variables (see 3.12), and B.) where g is 
the spherical radial coordinate r and the function F is an analytical function of a certain 
combination of the spherical angular variables (see 3.19). In both cases F and also g 
each solve Laplace's equation. Theses solutions are closely related to the possible forms of 
linearly polarized TEM waves.0 

Our method of solution is very particular and and we will employ it in further work. 
Finding the solution to equations (1.2), (1.4) is a first step in an attempt to solve three di- 
mensional hydrodynamical problems using potential methods. The procedure is as follows. 
The transformation from the three Cartesian variables variables to the functions g,F,F 
leads naturally to a metric function q^ 

ds 2 = dx\ + dx\ + dxj = q ij dy i dy j . (1.5) 

Where we have used y\ — g, y 2 = F, y 3 = F. First, the partial differential Eq. 1.4 
may be interpreted as constraints on the metric tensor: the appropriate relations among 
the elements of the metric tensor are found. Next with this metric one may calculate 
the Ricci tensor. In dimensions fewer than 4, the Riemannian curvature tensor and the 
Ricci tensor are linearly related one to another (see Landau & Lifshitz§, for example). 
Indeed, the vanishing of the Ricci tensor is both a necessary and sufficient condition that 
the coordinate system g, F, F be a diffeomorphism of the Cartesian x±, x 2 , £3 system. It is 
interesting that only in the very last step of the procedure, that one makes the correspon- 



dence of the g, F, F system with the Cartesian system. That is, the original vector pde 
is solved with respect to the natural coordinates of the problem, which are not in general 
the Cartesian variables. 

In section 2, we show how the solutions to 1.1 are related to solutions to Euler's 
equations. In section 3, we completely characterize the solution to 1.1 by the method 
outlined above. Since the method is novel, we show what the method might look like as 
applied in solving two-dimensional incompressible, potential flow problems in Appendix B. 
Of course, in the latter case, the solutions may be derived straightforwardly by other well 
known means; the discussion is included since it isolates the method particularily well (for 
some already familiar problems). The reader is strongly encouraged to read Appendix B, 
before beginning section 3. 

Section 4 discusses what sort of boundary conditions the complex Trkalian fields may 
satisfy. The type A fields do not satisfy natural boundary conditions, however the type B 
(singular) fields may have no normal component to a bounding surface. 

Section 5 discusses how the work presented here meshes with Bj0rgum0 and Bj0rgum 
and Godal's^l work on Beltrami fields. Section 6 is a discussion and conclusion. In ongoing 
work, we discuss solutions to other vector pde's: including linearly polarized transverse 
electric field§ (TE) and TEM@ solutions to the the electromagnetic wave equations. The 
main purpose of this article is to introduce a potential method capable of constructing truly 
three dimensional vector field solutions to pde of mathematical physics. 

A glossary of various terms is appended at the end of the manuscript. 
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2 Hydrodynamics and The Complex Trkalian Field 



Trkalian fields yield solutions to Euler's equations for an incompressible inviscid fluid. 
These latter equations may be written: 

where the density p = constant; and V • v = 0. A pressure function may be found (at least 
locally), if and only if the vorticity equation is satisfied: 

§ + Vx (fxtf) = 0, (2.2) 

— # — * 

where £ = V x u . Thus any time independent vector field, satisfying £ x v = is a solution 
to Euler's incompressible equations .0^0 We shall call these flows Beltrami flows (as is in 
keeping with Bj0rgum). A subset of these flows are the Trkalian fields, satisfying £ = Xv, 
where A is a constant, which under suitable scaling and inversion we can always take to be 
unity. 

Consider therefore the Trkalian flows which solve the differential equation, 

f=Vx{7={7. (2.3) 

Trkalian fields are divergenceless and obey a linear superposition principle. Each component 
of the vector field satisfies the Helmholtz eigenvalue equation, 

Avi + v i = 0. (2.4) 

In the literature therefore, Trkalian fields are often represented as a sum over Fourier modes, 
using wavevectors with unit norm: 

(ffj + ifexo^, (2.5) 

where k ■ = and k ■ k — 1. It is our contention that the Fourier type representation 
(2.5) is not necessarily the most useful representation of Trkalian fields. 



Time-independent Trkalian fields (2.3) clearly solve the vorticity Eq. (2.2) for the 
incompressible Euler's equations. Trkalian fields multiplied by the time decaying factor e~ vt 
are solutions to the incompressible Navier Stokes equations (Eq. 2.1 with a term uAv 
included on the RHS, where v represents the viscosity: see 

For this paper we are interested in those vector fields which solve the Trkalian equation 
2.3, and have a "topologically dual field," w, which also solves the Trkalian equation. The 
relations between vector fields and their topological duals have recently been investigated 
by one of the authorsjil some results are collected in Appendix A. In summary, a vector 
field, v, and a dual field, w, are related in such a way that 



for some g, F and F, where g, Fi, F 2 are all real functions. As pointed out in the Appendix 
any real vector field may be represented as the real part of a complex vector 
field, written as 2.6. That is c (as well as w) may be expressed with precisely the same 
three functions with which one may express v. Among other properties, a vector field is 
always orthogonal to the curl of a dual field, and a field and a dual field both have precisely 
the same value of the helicity H = v- Vxv = w- 'Vxw. A given positive helicity vector 
field does not in general have a unique dual field (although a suitable gauge fixing condition 
could probably be defined.) 

We call the simplest possible solution to (1.1) and (1.2) "the simple helical field": 



c = v + iw 



(2.6) 



v = (sin z, cos z, 0), 



(2.7) 



with 



w 



(— cos z, sin z, 0), 



(2.8) 



c = 



v + iw 



0) = e tz V(x + iy) ■ (-z). 



(2.9) 



6 



Eq 2.9 is obviously of the form of (1.1). Notice that v and w are orthogonal and both have 
unit helicity (implying c- c = where dot is the usual scalar product). For the Trkalian 
vector fields, we may very roughly think of a dual field as a rotation of each vector around 
the screw axis (in this case the z-axis), where the screw axis is along the Vg direction. 

3 Solution to the Complex Trkalian Equation 

The reader is urged to read Appendix B so that the logic used in our approach will 
be completely clear. We begin with the complex Trkalian equation: 

V x c = c. (3.1) 

Suppose in addition that 

V x c = iVg x c, (3.2) 

— * 

where g is some real function. Eq. 3.2 implies that c — e tg V F, for some complex valued F. 
That is, c is a zero-helicity field of a very special type. (In general, c — GVF, where F, G 

— * — * 

are complex, yields c • V x c — 0.) Substituting the relation, c — e tg VF, into 3.1 yields 

fgxVF = -iVF, (3.3) 
VgxVF = +iVF, (3.4) 

where F is the complex conjugate of F. We wish to find the components of the tensor, = 
^rf^r, where we define y± = g, y 2 = F, y 3 = F. First we note that this tensor and the 
tensor q lk = f^f^- are matrix inverses. Next from (3.3), for example, one immediately 
notes that VF • VF = 0. Thus q 22 = 0. Similarily we find g 33 = 0, g 13 = q 12 = 
and q 11 = 1. The last result comes from taking the inner product of (3.3) and (3.4) and 
simplifying. Finally define H' = Vy l ■ Vy 2 x Vy 3 , which turns out to be a pure complex 
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entity. As is always true, 


we have det q ab = 


i n ) 


Putting 


these together, we may write 






' 1 














q ab = 



_ 




iH' 


iH' 



J 


(3.5) 


with matrix inverse 




" 1 














q a b = 







1 

iH' 


i 

iH' 






(3.6) 



1. This is equivalent to 

ds 2 = dg 2 + 2e 2J dFdF, (3.7) 

if we define e 2J (as yet undetermined) to be The no n- vanishing components of the 

metric tensor have been identified. This completes the first step of the procedure. 

2. At this point we have just calculated the metric tensor. We must ensure however 
that the various coordinate relations encoded into the metric tensor are consistent with a 
flat (Euclidean) three-dimensional space. We begin with the line element 3.6 (or 3.7) and 
calculate the Ricci tensor. First 



pi = 1 im ( d( lmk dqrrd _ ggta 

kl ~2 q [dyi dy k dy m 



(3i 



Then 



be 



■pa 

1 bc, a 



ra I "pr "pa 
ba r ~T 1 hr L 



-pa pr 

br ca 



ba,c ' be ra 

' ' J\2 —Jl3 

-e 2J (J n + 2J 2 ) 




-2 (J 2 + J x 



2 J. 



23 



(3.9) 



where the subscripts following commas denote derivatives with respect to g, F and F re- 
spectively. The Ricci tensor R is symmetric, thus we have written down only the upper 
triangle in 3.9. Fortunately, there are symbolic manipulation programs available that allow 
one to calculate the Ricci tensor from the metric tensor in a matter of moments. 

The curvature tensor vanishes iff the g, F, F system is a diffeomorphism of the Carte- 
sian system. In 3 dimensions this condition is satisfied iff the Ricci tensor vanishes.il Thus 
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since Ji 



0, we need: 



J = In G(g) + In M(F,F). 



(3.10) 



Since Jf 



+ Jn = 0, then up to scaling and translation we have 



G — l/y/2 

R G = g/V2. 



or 



(3.11) 



Conditions A and B ensure that all the components of the Ricci tensor (3.11) vanish except 
for the (2,3) component. 

The case A leads to (lnM) 2 3 = (since J\ = in 3.11) implying 2J = — ln2 + 
In W'(F) + In W'(F), for some W, W which are arbitrary functions of F,F respectively 
(the prime ' denotes derivative). Then 2e 2J = W'(F)W'(F). Thus the line element may be 
written as ds 2 = dg 2 + dVTdW 7 . Thus we may define x and y by W = x + iy so that we 
identify g with the Cartesian variable z. That is, we may pick axes so that g = z and F is 
any complex analytic function of W = x + iy. Thus for case A . ds 2 = dz 2 + dx 2 + dy 2 , and 



For the simple helical field mentioned in section 2, we would have F(x + iy) = —i(x + iy) 
and c — e lz (—i, 1, 0). Curiously it seems to have been Beltrami0 who first wrote down that 
the Trkal equation was solved by the real part of 3.12; unfortunately this line of thinking 
has never been developed thoroughly until this article. 

The case B We wish the 23 component of the Ricci tensor to vanish: 



c = e iz V (F(x + iy)) = e iz F' ■ (1, i, 0). 



(3.12) 



e 2J (J u + 2J 2 ) + 2J 23 = 0. 



(3.13) 



We substitute J = -\ In 2 + In g + In M(F, F) (this ensures that all the other components 
of the Ricci tensor vanish) into 3.7 and 3.15 giving 




(3.14) 
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with 



, <9 2 lnM 

- M=4 » (3 - 15 > 



As one might expect from the line element 3.14, we can relate the case B to the study of 
the spherical coordinates. 

Equ. 3.15 has fascinating properties.@~0 It arises in the study of heat propagation, 
movement of plasmas and traces its origins back to Emden, Darboux,@ and LiouvillJl as 
well. In different contexts, it bears the names as Emden-Fowler, Arrhenius and Thomas- 
Fermi equations. First we note that it is equivalent to the two-dimensional equation = 
-e* (using M = -^e^ /2 ). Second we note that if M(F, F) is a solution to 3.15 then \A'(F)\M(A, A) 
is also a solution to 3.15 for any arbitrary complex analytic function, A(F) (or ip(F,F) is 
a solution implies that -> i/;(A(F), A(F)) + \a(A(F)) + \n(A(F)) also is a solution). 

Thus introduce A, A so that 

M 2 = A'(F)A'(F)/ cosh 2 (^-^). (3.16) 

This solves 3.15 for arbitrary complex analytic functions A, A. That any solution to 3.15 
may be expressed in the form of 3.16 was originally given by LiouvilldH (who expressed 
his solution with respect to a variable corresponding to the exponential of the function A.) 
For example, setting A = 6(lnc + lnF), where b, and c are real constants yields M = 
2b/(R((cR) b + (cR)~ b )), where R 2 = FF. Solutions of this type correspond to the much 
discussed radial solutions in flame propagation.il On the other hand setting A = —ib(\nc + 
InF), yields M — 1/(R cosh ( arg F)), where arg takes the argument of a complex number. 

Next define cos# = tanh and <fi = i(A — A)/2. Then it is not difficult to show 
that M 2 dFdF = d6 2 + sin 2 Odd) 2 . Thus 



ds 2 = dr 2 + r 2 (d6 2 + sin 2 6d(f) 2 ), (3.17) 
and we recover the spherical coordinate system. Since A was a completely arbitrary function 
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of F, then conversely F is a completely arbitrary function of A. So due to the definitions 
of 9, we have 

-Spherical = F(-arc tanh(cos6 l ) + 1(f)). (3.18) 

We should require also that either F be a 27r-periodic function of 0, or a (complex) scalar 
multiple of the identity function. This ensures that theresulting vector fields are single 
valued. Notice that F is anihilated by the operator <% + i-^gd^, which may be used to 
factor the angular momentum operator. And so 



c = 



9 + i$\F'. (3.19) 



rsin# 

A typical case is given in Figures 4 and 5. Choosing F(a) = cos(— m) and using 3.18 
and 3.19 yields v = rs 7 n a e |#[cosr cos0cos# + sinr sin0] + 0[cosrsin0 — sin r cos cos 9] j . 
Thus even a very simple function chosen for F may yield a rather complicated expression 
for v. 

In summary, there are precisely two cases: A. the Cartesian case (3.12), where the 
screw direction is in the direction of a Cartesian variable and F is an analytic function 
of the other two Cartesian variables, and B. the spherical case (3.19), where the screw 
direction is in the radial direction of a spherical coordinate system and F is an analytic 
function of the variable (—arc tanh(cos6 l ) +i<f>), with the additional requirement that the 
resulting vector field be single valued. 

It is straightforward to calculate the helicity for either cases A or B. We find 

Case A: h = v 2 = -VF ■ VF = F'F', 

2 

CaseB: h = v 2 = -VF ■ VF = F'F'/ (r 2 sin 2 9), (3.20) 

2 

where F' depends on either A. x + iy or B. —arc tanh(cos 9) + i<f>. In both cases AF = 0, 
the 3-dimensional Laplacian of F is zero for either class. In case A, actually the two- 
dimensional Laplacian is of course zero, and in case B, the "angular momentum part" of 
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the Laplace operator acting on F is zero. 

As a final remark, notice from 3.20 that there are no oscillations in the helicity (equiv- 
alently velocity magnitude) in the direction of increasing g, (that is z or r.) This is related 
to the fact that c is a complex null vector c • c = , so that only the cross terms v 2 = 
survive in calculating v 2 . 

From 3.20, we notice that nothing precludes that the helicity should vanish at excep- 
tional points. Note also that the integrated helicity over any infinite domain of R 3 (for the 
complex Trkalian flows) will always be infinite. In that sense, these solutions lack a certain 
natural global character. 

4 Boundary Conditions 

Bj0rgum has shown that in regions free from singularities, that if either v ■ n = 
or v x n = 0, where v is a solenoidal Beltrami field and n is the normal to some bounding 
surface, that it necessarily follows that the vector field v must vanish identically. We are 
then able to discuss which types of boundary conditions the Complex Trkalian Fields 

may obey. 

Since the type A complex Trkalian fields (with screw direction along a Cartesian axis), 
are solenoidal, they cannot satisfy the aforementionned natural boundary conditions. The 
type B Complex Trkalian fields always have a singularity interior to any region containing 
the origin so that the rigorous result of Bj0rgum does not apply. Moreover these fields 
always have a vanishing component in the radial direction, so that one may pick a boundary 
( r = a constant) that encloses a simply connected region in such a way that the vector field 
lies tangent to such a surface. However, there is a line singularity passing through the center 
of the sphere, leading to an infinite value of the integrated helicity. In conclusion then, 
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we see that no finite amplitude complex Trkalian field in a simply connected domain may 
satisfy the physically interesting boundary conditions tf-n = on some surface. In multiply 
connected regions, Trkalian fields satisfying physically interesting boundary conditions may 
exist. 



5 Bj0rgum and Godal's Work on Beltrami Fields 

In this section, we discuss Bj0rgum et. al.'s work (1948, 1951, 1952) and intepret his 
results in the same framework as ours. In his first long paper on Beltrami fields,0 Bj0rgum 
points writes that the Beltrami vector fields are precisely those vector fields which may be 
written as: 

Vgi x Vq 2 



ft 

where ft is the torsion function. This leads to 



<?iVg 2 + Vg 3 (5.1) 



df df 
dq a dq b 



R S 

S T ?1 § 

n a n n 

U qiH h 



(5.2) 



where R, S, T, M are themselves functions, and T may be found from: 



(The determinant of the matrix 5.2 is 1/H 2 , where H = Vqi ■ Vq 2 X Vg3 is the Jacobian 
of the transformation from the g's to the Cartesian variables.) One may demonstrate Eq. 
5.2, using the identity: 

dr = Vgi x Vg 2 

dq 3 Vqi ■ Vq 2 x Vg 3 ' 
which holds for any three functionally unrelated functions q±, q 2 , q 3 - Using 5.4, the Bel- 
trami equation (Eq. 5.1) may be written 

df „ ( df df df df\ ,„ „ x 

^^^ x ^ + ^ x d- (55) 
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Equation 5.2 then follows directly. Bj0rgum, however does not ensure that the above metric 
tensor given by 5.2 yields a zero curvature tensor. This must be the case if the q\, q 2 and q 3 
coordinates are to be themselves functions of Cartesian coordinates, as we argue in our 
paper. 



If one uses the representation given in Eq. A. 7: 



1 



v = -(e^Vy 2 + e-^Vy 3 ), 



(5.6) 



where y 1 = q 2 is real and y 2 and y 3 are complex conjugates, then Bj0rgum's statement is 
equivalent to the following. The Beltrami fields are precisely those vector fields where: 



df df 
dy a dy b 



where 



B 


-{e iyi A 


\e~ iyi A ' 




-ie iyi A 


e 2iyi D 


P 


(5.7) 




P 


e -2iyi D 


ab 


B = 1 


±A 2 ttH 
HttR ' 






n 

D = 


-RH 

AH ' 






P = a 


+ RH 
AH ' 




(5.8) 



are real functions. Eq. 5.7 may be used as a starting point for the Ricci tensor calculation. 
Given Q and H, there are two functions (say R and S) that characterize the matrix (5.2). 
Similarily, two functions characterize the matrix given by 5.7 (say R and A). The function A 
is related to S by 

A = S + Rq 3 . (5.9) 

The determinant of the matrix 5.7 is [D + P][B(D - P) + a 2 /2] = = This is 
because H' = Vyi • Vy 2 x Vy3 = — 2iVgi • Vq 2 x Vg3 = —2iH as may be found from 
calculating t?-Vxi; from 5.1 and 5.6. To proceed from 5.2 to 5.7, one needs 



df df dqi df df dqj 



df df 



dy a dy b dy a dq { % dy b dq t dqj 
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(5.10) 



with 



dy a 



1 -9i 



fe^ 1 ±e" iyi 



(5.11) 



The above matrix [Z] may be found straightforwardly by using yi = q2, 2/2 = (<l3 + 
iqi)e~ iq2 , ys = (qs — iqi)e +iq2 . The Beltrami equation in terms of the y variables then may 



be written 



e iyi i 



df df ( 

+ ie~ iyi — = n I e iyi 



dy. 



dyi 



df df 
dy 3 dyi 



+ e 



df df 
dyi dy 2 . 



(5.12) 



For complex Trkalian fields it is straightforward to show in that case that one must 
have A = D = and B = 1 (which specifies the torsion function). This matches the 
coefficients of e m and e~ tyi above in Eq. 5.12, or equivalently, discards the oscillating 
terms of the matrix, 5.7. At that point one is left with the matrix described below Eq. 3.6: 



df df 
dy a dy b 



1 





i 

iH' 

o W 



(5.13) 



ab 



The curvature condition then ensures a complete solution. In future work we will proceed 
to examine less restrictive cases than the complex Trkalian case and calculate curvature 
more directly from Eq 5.7. One might suspect that the complex Trkalian fields hold a very 
prominent place in the theory of the Beltrami fields, since they seem to take the elements of 
the matrix 5.7, which are most robust in some sense. As we show in later work the complex 
Trkalian fields may be shown to yield precisely the linearly polarized TEM solutions0 of 
electromagnetic wave theory. 



6 Conclusions - Discussion 



In this paper we have introduced a new method for constructing true three dimensional 
vector field solutions to vector pde's of mathematical physics. As opposed to axisymmetric 
flows, which are the dominant solutions to fluid mechanical equations listed in texts, true 
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three-dimensional flows should have a finite value of helicity v ■ V x v. It is clear that the 
Clebsch representation is then the most convenient: v = fVg + V/i. That is because non- 
zero helicity fields may always be expressed in such a manner^ and can never be expressed 
with fewer than three Clebsch functions. Given the representation, however, the question 
remains how to construct these functions (f,g,h). 

Our approach then is to represent the appropriate vector field with the appropriate 
Clebsch functions and then to construct the relations between the elements of a metric 
tensor, q^, using the original pde. After having constructed a metric tensor consistent with 
the equations one wishes to solve, one computes the Ricci tensor. If the Clebsch functions 
are to be functions of the Cartesian variables, then (in three dimensions) the Ricci tensor 
must necessarily vanish. After solving the partial differential equations for the vanishing of 
the Ricci tensor, one finally searches for the class of transformations of the metric tensor 
back to Cartesian coordinates. By this means, the integrity of the good coordinates - the 
Clebsch functions - remains intact. 

In particular in this manuscript, we have undertaken the first steps in a systematic 
classification of Beltram&@ and Trkalian fields which are thought to be important con- 
stituents of high Reynolds number flow fields. For example the so-called ABC (Trkalian) 
flow has been extensively studied^ particularily for its chaotic properties.0~@ Note that, 
although it is always possible to write down a formal solution to the Trkalian equation as 
in equ. 2.5 in terms of Fourier modes, the use of Clebsch functions, which we advocate 
here, has many advantages. For example, when one Fourier decomposes a problem one 
implicitly assumes that the Cartesian coordinates are a natural coordinate system for the 
problem. As we have discussed above this is not necessarily the case. 

In this particular article the analysis proceeds by considering the Beltrami fields as 
real parts of complex vector fields, and then considering complex vector fields of a cer- 
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tain type. Complex vector fields play a prominent role in many areas of mathematical 
physics including electromagnetism (see Bateman@) and minimal surfaces (see Struikil or 
Osserman0) In this paper, we examine complex vector fields, which satisfy the Trkalian 
equation (1.3), and may be written in the form e ?9 V-F, with g real and F complex (that 
is, we restrict the problem enough in order to give a complete solution). Both the real and 
imaginary parts of such a vector field will then be Trkalian. The two classes of solutions to 
the complex Trkalian equation with the complex vector field of the above form are given 
in section 3. The class of solutions where v is a 2-d potential field with a screw direction in 
a Cartesian direction has been studied before (Bj0rgum's uniplanar Trkalian flows). The 
class of solutions where v is a "zero angular momentum field" with a screw direction in 
the spherical radial direction does not appear in the Beltrami literature. Examples of both 
types of fields are presented in the figures. 

Although the Clebsch representation of the vector field has found fruitful application in 
fluid mechanicsii - ^, thermodynamics and more generally in variational principles the 
authors still maintain that this representation has been underexploited. We demonstrate in 
this article that the method we have used employing Clebsch functions has much practical 
merit in producing new solutions to partial differential equations. The same method will 
be used in further articles on fluid mechanics, and wave equations. 

The author would like to express his deep gratitude to R. M. Kiehn at the University of 
Houston for teaching him about differential forms, the intrinsic merit of Clebsch functions, 
and a topological approach to Physics and for sharing his results prior to publication.0 The 
author also thanks Professor Giles Auchmuty of the University of Houston for discussions 
and for pointing out several references. 
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A Vector Fields and Their Topological Duals 



All the following arguments are local for a region where the helicity is of constant sign 
(non-zero). 

Proposition 1 Every vector field, v, may be written as the real part of a complex vector 
field of the following form: 

c = e !S VF, (A.l) 

where g is a function with domain 1Z 3 and range 1Z, and F has domain 7Z 3 and range the 
complex numbers. 

Proof: One knows that it is always possible to construct the Monge potentials /, g, h so 
that a vector field v may be expressed in the decomposition due to ClebschS 

v = ffg + Vh. (A.2) 

Now define 

F = (h + if)e- i9 , (A.3) 

F = (h- if)e +i9 . ( A.4 ) 



Or equivalently, 



Now we may write v as 



from which 1 follows. 



/ = -^(Fe^-Fe^), ( A.5 

h = X -{Fe i9 + Fe~ i9 ). ( A.6 



v = ^(e i9 VF + e~ i9 VF), ( A.7 
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Proposition 2 The real part of c, v, and the imaginary part of c, w satisfy: 

V x v + Vg x w = 0, ( A.8 ) 

VxiH-VjxtJ = 0, ( A.9 ) 

vxw^O (A.10) 

We say that any two vector fields that satisfy A.8 - A. 10 for some function, g, are said to 
be dual to one another with respect to g. 

Equs. A.8 - A. 10 are easy to show beginning from 

w = -hVg + Vf. ( A.ll ) 

In particular it may be shown from the above that dual vector fields satisfy 

w-(Vxt)) = 0, (A.12) 
v-(Vxw) = 0, (A.13) 
v-(Vxv) = w-(Vxw) (A.14) 

V x ( % X Z\-^ X :\ v- %^P^- w - ^— vxw) = 0. ( A.15 ) 

\(v X w) ■ (v X W) [V x w) ■ [V x w) [V x w) ■ [V x w) J 

— * 

The terms in the parenthesis of A.15 constitute a solution for —Vg from A.8 and A.9. 

A complete classification of vector fields in this fashion does not seem to be known. 
For example, one can also take a vector field in TZ 4 and write it as the real part of a zero 
helicity vector field with domain the complex numbers. One can also write a vector field 
in 1Z 8 as the real part of a zero helicity vector field with domain the quaternions. For 
the dimensions 5, 6 and 7 however no obvious standard form exists, although clearly these 
vector fields may be thought of as the real part of zero helicity quaternionic fields of a 
restricted type (the same way that any real vector field may be thought of as the real part 

— * — * 

of a complex vector field satisfying, V x c = iVg x c.) 
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Dual concepts are prevalent in mathematical litaerature. Typically one hopes to rep- 
resent a function (or vector of functions) as the real part of a complex function with slightly 
nicer properties. The complex function uses the same functions as the real part as build- 
ing blocks. Thus one writes a harmonic function as the real part of an analytic function 
(which depends on a single variable, but where this variable belongs to the complex field). 
Similarily one writes the panharmonic function as the real part of a fi regular function (see 
Duffini). 

B The Proposed Method Employed in a Familiar 
Problem 



We give a very simple example of the general procedure we are outlining. Consider 
the case of incompressible, time independent, irrotational flow past an obstacle in two 
dimensions. We wish to study a potential function, 0, a stream function ip, and a complex 
potential W = <fi + iip. Let's try to write a differential element of length in terms of the 
functions and if). Since <fi and ip are respectively the real and imaginary part of the 
complex potential W, they satisfy the Cauchy-Riemann relations: |^ = — t^; = t^, or 

x- J ax ay 1 ay ox 1 



more compactly 



V0 = z x Vip. 



B.l 



Now one may straightforwardly calculate 



V0- V0 
V0- 



0. 



B.2 
B.3 



For shorthand, let y 1 



, 2 = ^. Then g lk = ^M 



9ij 



dx p dx p 
dy l dyi 



a 2 
a 2 



1/a 2 
1/a 2 



. The matrix 



B.4 
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is the matrix inverse of g lk . 

1. The first step of the calculational procedure is then complete: for the equations of ir- 
rotational, incompressible flow (Cauchy-Riemann conditions), the appropriate non- vanishing 
components of the metric tensor are identified. We are then left with a line element (or 
metric tensor) given by 

ds 2 = a 2 (d0 2 + d^ 2 ) , ( B.5 ) 

where a = a(0, ■?/>). Of course, B.l immediately implies B.5 due to conformality of a 
complex anayltic transformation; one needs virtually no calculation. In the approach taken 
here, flow equations like B.l are interpreted as geometric constraints (on the 
metric function). Note that in his papers Bj0rgum calculated correctly the non- vanishing 
components of the metric tensor. He never ensured however that the curvature was zero. 
Applying this spirit to the current problem, Bj0rgum would have halted at this point. 



2. Now we calculate the Ricci tensor given the metric function B.4: 

" a 2 



9a 

Let A = In a. Then 



a 



2 



B.6 ) 



r i _ 1 im ( d 9mk dg mt dg M \ 

= A,i5 ik + A, k 5 u - A,i5ik , ( B.8 ) 

where the indices following commas denote derivatives. Now the Ricci tensor is given by 

Rbc = IT - T£ + ^IJ^ra - r br r L = -(-4)11 +-4,22 



1 

1 



( B.9 

be 



Now we may express the functions <fi and ip in terms of the Cartesian variables if and only 
if the Ricci tensor vanishes identically. Thus we need the function In a to be harmonic with 
respect to the functions and tp. Actually this is quite obvious since 

ds 2 = dzdz = — - — dWdW, ( B.10 ) 

dW dW 
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where z = x + iy, so that d ^ dw In (^jf?^) — 0- The above "A" corresponds to In 

3. Boundary Conditions 

The actual choice for In a must be determined by the boundary conditions of an actual two 
dimensional flow problem. For flow past an obstacle, we wish to map (part of) the real axis 
of the variable W = <fi + iip to the boundary of the obstacle. The upper half plane is then 
mapped to the exterior of the obstacle in the upper half plane. Once a((j>, ip) is determined, 
then we may try to calculate the Cartesian variables by means of the eikonal like equations 

There is no new result outlined in the approach above, we merely demonstrate the 
method, which we use to construct Clebsch functions in this paper, to a more familiar 
problem. 

4. Examples 

Symmetrical flow past a unit disc with unit speed far downstream. We need a function 
mapping (part of) the real line to the boundary of a unit disc. This is well known: z = 



x + iy = e iarccos f = (W + VW 2 -A)/2 We then calculate a 2 



dz 



dW 

So ds 2 = a 2 (d0 2 + dip 2 ), where 0, ip are the real and imaginary parts of W . 



1+ w 



VW 2 -4: 



For symmetrical flow past a plate of width 4 aligned perpendicular to the flow along 



the x direction, we have z = y/W 2 — 4, and we calculate a 2 



dz 



dW 

For both plate and disc, there is a singularity in a when cf) — 2; ip — 0. 



v /(^2 + ^2_ 4 )2 + 16 ^2 ' 



Conversely given the function a 2 , (where In A is harmonic in if),<f>) then one may use 
B.ll to reconstruct the boundary. 

Also note that since the streamlines correspond to contours of constant ip, then roughly 
the function ip plays the role of an energy for a 1 degree of freedom Hamiltonian system, 
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whereas the function <p plays the role of a phase (or time) variable. 

As is well known, the harmonic functions are precisely those functions which serve as 
nice coordinates of orthogonal coordinate systems in 2-dimensions. They are the "isother- 
mal parameters''^ of the plane. In a similar sense, the (Clebsch) functions we solve for in 
this paper are also the good coordinates for the problem at hand: the building blocks of 
the complex Trkalian flows. 
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FIGURE CAPTIONS 



Fig. 1 A picture of the simplest Beltrami field (a uniplanar Trkalian field of the first kind). 
Here g — z and F = x + iy. The patterns are reminiscent of the cholesteric mesophase 
of liquid crystals. 

Fig. 2 Here are two planes of the flow g = z and F = cos(x + iy). Each arrow from the 
lower plane to the upper plane is rotated the same amount around the screw axis 
(the z-axis). In each plane, the vector field is potential flow (see Fig. 3) 

Fig. 3 The bottom plane of Figure 2. 

Fig. 4 A total of 714 vectors are plotted on the surface of the sphere showing the complex 
Trkalian flows for g = r and F = cos(0 + % arc tanh(cos 9)). Here r is fixed to be 0.5. 
The vector magnitudes are not drawn to scale. The true vector magnitudes are the 
fourth power of those represented in the figure. This is to keep the arrows at the 
poles from dominating the figure (where the vectors become singularily large). 

Fig. 5 A better look at one of the singularities of the complex Trkalian flow of Fig. 4. 
There are a total of 4 singularities. Two singularities lie at the poles (on the z-axis) 
and two more singularities lie on opposite sides of the sphere (at the same latitude). 
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TABLE CAPTIONS 



Table 1 Some elementary fields: Note the word complex has a different connotation in 
its usage in complex lamellar and complex Laplacian (where we mean "rescaled" or 
"conformally equivalent to) than it does in its usage in complex Trkalian. Complex 
Trkalian describes a vector field whose imaginary and real parts are both Trkalian. 
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